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Abstract 

We introduce a finite dimensional matrix model approximation to the algebra of 
^< functions on a disc based on noncommutative geometry. The algebra is a subalgebra 

of the one characterizing the noncommutative plane with a * product and depends 
on two parameters N and 9. It is composed of functions which decay exponentially 
outside a disc. In the limit in which the size of the matrices goes to infinity and the 
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noncommutativity parameter goes to zero the disc becomes sharper. We introduce 
a Laplacian defined on the whole algebra and calculate its eigenvalues. We also 
calculate the two-points correlation function for a free massless theory (Green's 
function) . In both cases the agreement with the exact result on the disc is very good 
already for relatively small matrices. This opens up the possibility for the study 
of field theories on the disc with nonperturbative methods. The model contains 
edge states, a fact studied in a similar matrix model independently introduced by 
Balachandran, Gupta and Kiirkguoglu. 



1 Introduction 



Noncommutative Geometry PJ El El Ej provides the possibility to describe also ordinary 
geometries approximating the fields defined on them with finite (matrix) algebras. The 
archetype of these fuzzy spaces is the fuzzy sphere Along the same lines other fuzzy 
spaces have been built [HI |H|. There is also the fascinating possibility that the very 
structure of spacetime at ultrashort length scales might be fuzzy, with a dramatic drop of 
degrees of freedom at short distances with its consequences on the ultraviolet structure. 
As additional example matrix models in noncommutative geometry also provide a way ^T] 
to ask questions (for the moment in a toy model) like the number of points or the average 
dimension of the universe, questions which was unthinkable to ask in the usual context. 

The fuzzy approximation is in some sense similar to the usual lattice one, in that a 
continuous action is substituted by a model with a finite number of degrees of freedom, 
thus enabling the possibilities to solve field theories with a approximate path integral 
formalism. In the limit in which matrices become large the approximation improves. 
However there are some fundamental differences. One advantage of the fuzzy approxima- 
tion over the usual lattice is that the basic symmetries of the space at hand are preserved, 
thus fuzzy spheres are SO (3) invariant, fuzzy tori are translationally invariant etc. The 
original spaces are recovered when the size of the matrices goes to infinity. In this sense 
noncommutative geometry, and fuzzy spaces in particular, naturally lead to matrix mod- 
els jSlini!, which can approximate both commutative and noncommutative spaces. The 
study of field theory on these spaces with numerical methods has started O EH Ej 
and seems very promising. In particular the authors of ^3] solve exactly, on the Moyal 
plane, a matrix model which has similarities with the one presented in this paper for 
the disc. The finite matrix model can also be useful as a cutoff for renormalization, as 
advocated for example in ^B]. The connections between finite matrix algebras coming 
from a noncommutative plane and finite geometries has been noted in J7| in the context 
of noncommutative Chern-Simons theory ^Hj as a "quantum Hall droplet", while the 
connection between the projectors we use and the disc appears in 0120]. A model which 
has several points in common with ours has been independently developed by Balachan- 
dran, Gupta and Kurkgiioglu [21] in a matrix model which represents the edge currents of 
a Chern-Simons theory on an infinite strip. The difference between our model and theirs 
lies mainly in the way the limit is performed. 

In this paper we present a fuzzy approximation to the two-dimensional disc. We has- 
ten to add that although the root of the approximation lies in noncommutative geometry 
(or even earlier, in quantization) the space of functions being approximated is a normal 
commutative space of functions. In the limit in which matrices become infinite the non- 
commutativity disappears 1 . Another important aspect of the model is that, while in the 

^dn case the field theory under consideration is a nonabelian theory there can of course be other sources 
of noncommutativity which will remain. It is also possible (depending on the limit taken) to describe 
theories on noncommutative spaces. 
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lattice approximation the functions, and hence the correlations, can be calculated only at 
a finite set of points, in the fuzzy approximation presented here approximate functions 
have values at all points on the plane. We will argue that the matrix model we present 
corresponds to a theory on a disc in a variety of ways. We start in section |21 with a brief 
reminder of the noncommutative plane in the form which is useful for our purposes, which 
makes evident the map between functions on the plane and infinite matrices. In section |3] 
we introduce the fuzzy disc, a particular finite matrix subalgebra and we argue that the 
corresponding functions have mainly support on a disc and decrease exponentially outside 
it. As the size of matrices increase the boundary becomes sharper. We also discuss the 
form that the ultraviolet-infrared mixing takes in this case and argue the presence of 
edge states. In order to show that the algebra we are proposing actually corresponds to 
a disc, in section 0] we introduce derivatives and Laplacians and compare the spectrum 
of the latter in the continuum and in the finite approximation. In section |S] we sketch 
an approach to field theories on a disc based on the fuzzy approximation. Although the 
construction is well suited for Montecarlo techniques, the correlation of a free massless 
scalar field can be solved analytically giving the Green's function in terms of the inverse 
of a matrix. We calculate it numerically and compare it with the exact case. The final 
section contains conclusions and outlook. Some technical details of calculations are in the 
appendix. 

2 The Noncommutative Plane as a Matrix Algebra 

We start from functions on K 2 , with coordinates x and y. Consider this plane to be like 
a 'quantum phase space', that is that we "quantize" x and y and give them a nontrivial 
commutator: 



It is convenient to consider the plane as a complex space with z = x + iy. The quantized 
versions of z and z are the usual annihilation and creation operators, a = x + iy and 
= x — iy with a slightly unusual normalization, so that their commutation rule is 



The parameter 6, which has the dimensions of a square length, plays the role of h, but it 
has no physical meaning, much like the distance between sites in a lattice approximation. 

This "quantization" associates operators to functions. There are however ambiguities 
in this association and in order to give a unambiguous procedure we introduce a map 
which associates an operator on an Hilbert space to a function on the plane. This is a 
quantization procedure and it was essentially introduced by Weyl long ago [22]. Given 
a function on the plane we define a map fie which to the function ip(z, z) associates the 
operator <p defined as follows. Given the function ip(z, z) consider its Taylor expansion: 




(2.1) 



[a, a'] — 9 . 





(2.3) 
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to this function we associate the operator 



n 9 (<p):=<p = <pS* a r n a n , ( 2 - 4 ) 

thus we have "quantized" the plane using a normal ordering prescription. To be precise 
this map is not the one introduced by Weyl (which associates hermitian operators to real 
functions) but can be cast in the same language since 

fyfoO = / ^(u, u) e -V , (2.5) 

where <p(u,u) is the Fourier transform of (p. 

Because of the similarity of its integral expression to the standard Weyl map we will 
also refer to it as the weighted Weyl map. The map Qg is invertible. Its inverse is a 
variation of the Wigner map (23] (the inverse of the Weyl map). It can be efficiently 
expressed defining the coherent states: 

a\z) = z \z) , (2.6) 

then it results 

^e l i.v) = <f(z,z) = (z\ (p \z) . (2.7) 

There is another useful basis on which it is possible to represent the operators. Consider 
the number operator 

N = a)a , (2.8) 
and its eigenvectors which we indicate 2 by |n): 

N|n) = n6\n) . (2.9) 

We can then express the operators within a density matrix notation: 

oo 

= l m ) ( n \ ■ ( 2 - 10 ) 

m,n=0 

The elements of the density matrix basis have a very simple multiplication rule: 

\m) (n\p) (q\ = 5 np \m) (q\ . (2.11) 
The connection between the expansions ()2.4j) and ()2.10j) is given by: 



^y/(n + l)9\n) (n + l\ 



n=0 
oo 



a 

n=0 



' -= ^2y( n + 1 ) d \ n + 1 ) ( n \ ■ ( 2 - 12 ) 



2 There is a possible notational confusion here. While it is true that the coherent state for z — is the 
same as the eigenvector of N with zero eigenvalues, it is not true that the coherent state \z) \ z=n is the 
same as \n). Since we never consider coherent states with integer values we refrain from the introduction 
of another symbol for the eigenvectors of N. 



Some useful conversion formulas are collected in appendix El Note that with our conven- 
tions x, y, z, a and a* have the dimension of a length, N and 6 those of a length square, 
while density matrices are dimensionless. Applying the dequantization map ()2.7|) to the 
operator (p in the number basis we obtain for the function (p a new expression, analogous 
to the Taylor expansion ()2.3|) in terms of the coefficient (p nm , that is 

I 2 00 z m z n 

^"- VnM^ 1 (2 ' 13) 

where we used the expression (|A.3j) . The ordinary Wigner map of these density matrix 
elements (the Wigner functions) can be found in [25] . 

The maps Q and yield a procedure of going back and forth from functions to 
operators. Moreover, the product of operators being noncommutative, a noncommutative 
* product between functions is implicitly defined as 

(<p*<p')(z,z) = n- 1 (n(<p) n(<p')) . (2.14) 

This product (which is a variation of the Moyal-Gronewold product [23] at the basis of 
deformation quantization) was first introduced by Voros [2H] • Its differential expression is 

(if * if') (z, z) = e ed *' d *"(p(z')<p(z")\ z=zi=zn = tp<p' + 6d- z <p'd z <p' + 0(6 2 ) . (2.15) 

We will indicate the algebra of functions on the plane with this product as Ag. 

In the density matrix basis, because of (|2.11|) . the product (|2.15|) simplifies to an 
infinite row by column matrix multiplication: 

oo 

(f * V') m n = PmkV'kn ■ ( 2 - 16 ) 
k=l 

Although the map Q is not defined for 8 = 0, using equation ([2. 15)1 . it is easy to see 
that, when 6* — > 0, the * product goes to the ordinary commutative product. A word 
of caution, the various maps and products defined here have domains and ranges which 
are not identical. While the standard Weyl map maps Schwarzian functions into Hilbert 
Schmidt operators, for the weighted Weyl map ()2.4|) this is not always the case. 

Using relation (|A.3|) or directly ()2.13|) it is easy to see that also: 

„ CO 

/ <Pz<p(z,z) =7T0Tr$ = 7r0^ nn , (2.17) 

n=0 

where we have introduced the matrix $ with components (p nm . 



3 The Fuzzy Disc Subalgebra 

The main point of the previous section is that the map Q and its inverse provide a manner 
to associate to each function, ip an infinite dimensional matrix $. The price to pay is 
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that the commutative product of functions gets deformed through a parameter 9. In this 
section we define the subalgebras (with respect to the Voros * product), of finite NxN 
matrices. Considering the functions whose expansion (j2.13|) terminates when either n or 
m is larger than a given integer N, it is immediate to see that these functions are closed 
under * multiplication. In the limit of N — * oo, with 9 fixed, the noncommutative plane 
is recovered. 

The subalgebra can be obtained easily from the full algebra of functions via a projection 
using the idempotent function introduced in a similar context by Pinzul and Stern [2U] : 

N N 2n 2 

Pe N = E (-1 n) (n| z) = £ — -e~^ , (3.1) 

n=0 n=0 U -° 



where we have used the polar decomposition of z = re 1 ^ and (|A.3[) . The finite sum may 
be performed thus yielding 



jV _ r(jv + i,y) 
r ° ~ r(Ar + 1) ' 

where T(n, x) denotes the incomplete gamma function. In this notation it is then clear 
that, in the limit N —>■ oo and 9 — > with 

i? 2 = iV^ (3.3) 

fixed, the sum converges to 1 if r 2 /9 < N (namely r < R), and converges to otherwise. 
It has cylindrical symmetry since <fi does not appear. For N finite the function vanishes 
exponentially for r larger than R (see figure Q). 




Figure 1: The function for N = 10 2 . 



Already for = 10 3 it is well approximated (see figure EI) by a step function. 

The function Pg is a projector of the algebra of functions on the plane with the * 
product: 

P? * P» = P e N , (3.4) 



Approximate Identity 
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Figure 2: Profile of the spherically symmetric function P g for the choice R 2 
and N = 10, 10 2 , 10 3 . As N increases the step becomes sharper. 



N9 = 1 



and the subalgebra is defined as 

X = P s N *A e *P e N . (3.5) 

Cutting at a finite N the expansion provides an infrared cutoff. The cutoff is fuzzy in 
the sense that functions in the subalgebra are still defined outside the cutoff, but are 
exponentially damped. For example a Gaussian of width 9 (corresponding to the element 
|0) (0|) when projected, is mapped into itself, and is nonzero on the whole plane. The map 
for a normalized Gaussian centred around the origin but with width a > 9 also shows 
the infrared cutoff. Let us investigate the effect of P g N on some functions. We use the 
notation 

IT» =PF*<p*lf= £ ywr4- Z Z (3.6) 

m,n=l \/m\n\9 m+n 

and talk of "fuzzy functions" for the projected functions so that P g N = 11$ (1) is the fuzzy 
identity. At the level of operators Pq corresponds to the projector 

Pe = E |n> ("I • (3.7) 

n=0 

We call fuzzy disc the space corresponding to the algebra A g ', which is isomorphic to the 
algebra of N x N matrices. This, as a matrix algebra, is the same as the one for the fuzzy 
torus, the fuzzy disc and any fuzzy spaces in general. In [2S| the fuzzy sphere is indeed 
obtained with the use a projector on the Voros algebra of 1R 3 . What makes it a disc are 
the way to take the correlated limit of 9 and N keeping the dimensionful quantity R fixed, 
which we now pass to describe, and the extra structures we will discuss in section 0] 

The limit discussed here must be understood in an heuristic sense, since the algebra 
of functions on the disc is not approximatively finite one would have to consider it in a 
weak sense, like the way it is done in j2H]. A more rigorous analysis on the line of what 
is done in (201 would also be desirable. It must also be noted that different ways to take 
the limit give different commutative or noncommutative spaces, as discussed for example 
in (231111. 
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We first calculate the effect of the projector on a rotationally symmetric Gaussian 
centred at the origin of width a. Any cylindrically symmetric function tp(r) has a simple 
expansion: 

^(r) =e-Vj]^— , (3i 



n=0 



where </2 n can be calculated from (|A.4J) . In particular for the normalized Gaussian 

1 r2 



we have 



(p(r) = — e a 



1 / ^ 

= 1 



(3.9) 



(3.10) 



and the series can be summed to give 



2 JL r 2n 



I$(<p{r)) = e r o Y,^n-^ = e 



n=0 



_ £ r(A + 1,^(1-1 

7raT(N+l) 



(3.11) 



The map IT^ is just an infrared cutoff which in the limit becomes sharper and sharper. 
Equation (|3.11|) is plotted in figure El for 6 < a. 





Figure 3: Profile of the spherically symmetric function IT, 
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-e a 



for the choice R 2 



N9 = 1, N = 10 3 and for various choices of a compared with the unprotected function. 
Both functions are plotted, although inside the disc they are practically indistinguishable. 
The unprojected function is always the larger one. 



Things are very different if we try to localize the function at a distance smaller than 
y/9. Figure 0] shows what happens in such case: for the larger value, a = .66, the 
projected function is undistinguishable from the exact function (both are actually plotted 
in the figure). For the middle value, a = .56, with our numerical approximations, a small 
"bump" at r = R appears. Already for a = A96 the 'bump' has become a large Gaussian 
sitting at the infrared cutoff; the part close to the origin is still there, but it is quickly 
dwarfed by the infrared bump which grows very fast. For a ~ .4 it is already of the order 
of 10 17 . Keeping a fixed and increasing N, with R fixed, the bump disappears. The reason 
for this behaviour is that the factors (p n in (J3.1U|) become negative and smaller than -1. 
Therefore the individual terms of the series become larger and larger. Keeping the whole 
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series provides cancellations which do not happen truncating the series at a finite value. 
In fact the bump becomes negative for iV odd. This is a reflection of the fact that, with 
the weighted Weyl map, the operator corresponding to the Schwarzian function ()3.9|) with 
a < 6/2 is not compact, and hence it cannot be approximated by finite rank matrices. 




Figure 4: Profile of the spherically symmetric function H$ ( « j ) for the choice 
R 2 = N9 = 1, N = 10 2 and for various choices of a. 



In general the function ip is close to IT^ (<p) if it is mostly supported on the disc of 
radius R (otherwise it is simply exponentially cut) and if it does not have oscillations of 
wavelength smaller than y/9. It this case the function becomes very large on the boundary. 
This is a compact example of the ultraviolet-infrared mixing [31] which is one of the most 
interesting characteristics of noncommutative geometry. If we try to localize too much 
the function, unavoidably an infrared divergence on the boundary of the disc appears. 

There are however functions which are localized sharply near the edge of the disc. 
These are the edge states [H2j which play an important role in Chern-Simons theory 
and have been introduced in these matrix models by A. P. Balachandran, K, Gupta and 
S. Kiirkgiioglu in |2T] . and the discussion presented here is inspired by this work. The 
edge states are simply given by the highest one-dimensional projector: 

V edge = - d (z\N)(N\z)=e-ir w ^^ . (3.12) 

They are plotted in figure El for N = 10 and N = 100. It is evident that as iV increases 
they become sharper. 



4 Fuzzy Derivatives and Fuzzy Laplacians 

So far we have defined a projection from functions on the plane to a finite dimensional 
algebra of NxN matrices and discovered that, with an appropriate choice of N and 6 
we obtain a good approximation of functions supported on a disc. But, in the spirit of 
noncommutative geometry, if we want to talk of geometry, we have to define a Dirac 
operator or a Laplacian which give the metric properties of the system 0. To define 
derivatives thus we define a pair of operators which act on the same space as the elements 
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Figure 5: The edge states </> ed S e for N = 10 and N = 100. 



of A$ . Their commutators with a function define what we usually call the derivatives 
of a function. Of course we require that in the limit these two derivations converge to 
the usual derivatives of the disc. Combining them with 7 matrices would give the Dirac 
operator, to avoid discussing fermions we will concentrate on the Laplacian. 

The starting point to define the matrix equivalent of the derivations is: 

d z if = -(z\[a\tt(tp)]\z) 

9bV = jj(z\[aM<P)\\*) ■ (4- 1 ) 

The above expression is exact. Acting on an element of the subalgebra A$ the derivative 
takes the functions out of the algebra, a phenomenon not uncommon in noncommutative 
geometry. However 

d z (Ae)cA$ +1 (4.2) 

analogously for d s , so that the derivatives of functions of the subalgebra can still be 
considered finite matrices. Notice that 

0,<(¥O ^ itf&vO . ( 4 - 3 ) 

the equality obviously holding in the limit. We will use d z Il$ (<p) to define the Laplacian 
below. 

The fact that we keep a and a) operators on the full space (hence they are still 
infinite matrices) is crucial for the identification of the algebra of N x N matrices with 
the approximation to the disc. If we were just to truncate the matrices a and a\ their 
commutator would not be proportional to the identity. For the fuzzy torus the derivations 
cannot be expressed as commutators, while for the fuzzy sphere there are three derivations 
which are not however independent, but are connected by a Casimir. This goes just to 
say that the same algebra, with different structures, can approximate different spaces. 

Rotations are well defined, in fact the generator of rotations on the fuzzy disc is 
nothing but the number operator N introduced in (|2.8J) which commutes with Pg just 







as the generator of rotations on the ordinary disc, the angular momentum operator <9^, 
commutes with Pq . Just as the fuzzy sphere maintains the invariance group of the sphere, 
the fuzzy disc retains the fundamental symmetry of the disc. 

Note that the eigenvalue equations 



Xip 



(4.4) 



have no solution in the space of N x N matrices, just as in the commutative case trans- 
lations on the disc have no eigenvectors. Nevertheless the fuzzy Laplacian 
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can be defined. In particular consider the following matrix model: 



5 = ^/^^*( V V)=Tr$L(v 2 ) mnM ^, 



(4.5) 



(4.6) 



where (V 2 ) mriJ5(? is implictly defined by ()4.5|) . It is an operator acting on a finite space of 
dimension (N + l) 2 , and its eigenvalues can be calculated and compared with the exact 
commutative case. 

As is well known the eigenfunctions of the Laplacian on a disc with Dirichlet bound- 
ary conditions are the Bessel and the Neumann functions of integer order multiplied by 
exp(im0). The eigenvalues, E mn , are related to the zeroes of these functions. More 
precisely, we have 



E r , 



.zeroes. 



R 



(4.7) 



In figure |U] we show a comparison between the eigenvalues for the exact and approximate 
Laplacians for three values of N. The agreement is fairly good, a fuzzy drum sounds 



Eigenvalues of the Laplacian 
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Figure 6: The first eigenvalues of the Laplacian on the disc (dots) and the fuzzy Laplacian 
(solid lines) for N = 5, 10, 15. The lines corresponding to the three cases can be distin- 
guished because the agreement with the exact case improves as N grows. In the figure on 
the right the curve which interrupts is the one corresponding to N = 5, for which there 
are only 36 eigenvalues. 
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pretty much like a regular drum. Discrepancies start occurring for the 4N th eigenvalues, 
this is to be expected because AN/R 2 = 4/9 is the energy cutoff of the theory 3 . Most 
of the eigenvalues (but not all) of the fuzzy Laplacian are doubly degenerate, but the 
unmatched ones become sparser as N increases. 



5 Free Field Theory on the Fuzzy Disc: Green's func- 
tions 

The formalism we have developed lends itself readily for matrix approximations to field 
theories on a disc. For the real scalar case described by the action 

S = J d 2 z ipV 2 cp + + V(<p) , (5.1) 

we consider the fuzzy action 

Si = - Tr $V 2 $ + — $ 2 + V(<$>) . (5.2) 

7T 2 

We stress that this is an action entirely of finite dimensional matrices, which can be 
approached numerically, using for example Montecarlo techniques, a method currently in 
use for other fuzzy spaces [T2JE3E1- Here we will content ourselves in calculating the two 
points Green function for the free massless scalar theory. In this case the path integral 
may be performed yielding just the inverse of the Laplacian with Dirichlet boundary 
conditions which we have calculated in section 0] 

(<p(z, z)<p(* t z')) = G(z, z') = (z\ (V 2 )- 1 \z') . (5.3) 

In the continuous theory on the disc this can be calculated exactly using standard tech- 
niques of classical electrodynamics(see for example [33]) yielding 

G(z, z') = — In , , [ f ~ Z '[ - t . (5.4) 
v ' ; 2tt \z'\z\ - z\z\~ l \ y ' 

The fuzzy approximation is 



v ' - -^ ^{V- 2 ) mn pq zrz*z> m z' n 



U 9 {Z,Z)-4 2^ / .LULIflm+n ^ " ^-3) 



This expression can be evaluated numerically. In figures [7| and |S] we show a comparison 
between the exact Green function on the disc (|5.4|) and (|5.5jl . 

The agreement is, in our opinion, quite remarkable already for a limited number of 
points, the logarithmic divergence ha been smoothened by the ultraviolet cutoff, but apart 
from that the two functions are quite similar. The choice of different values of z' gives 
similar pictures. 

3 We use units for which Ti = c = 1 

1 1 



Figure 7: Comparison between the Green's function G(z, z') as function of z for z' = i/2 
On the left the exact function (with the singularity at z = z' truncated), on the right the 
approximated function for N = 20 . 
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Figure 8: The density plot of the Green's functions with the same parameters as in figure^ 
The exact function is on the left. 

6 Conclusions and Outlook 

We have seen as the fuzzy approximation to the disc works fairly well for the calculations 
we have attempted, which were based on the structure of the algebra and on the Laplacian. 
The first conclusion that can certainly be drawn is that the fuzzy disc, as a drum, sounds 
very much like a commutative disc even with an approximation based on relatively small 
matrices. We performed only one exact analytical calculation as a comparison with the 
continuum case, but the method can be extended to more involved calculations where it 
can be a very useful tool. Although the approximation is finite the resulting (fuzzyfied) 
functions, correlation functions etc. are defined for all points on the disc, and rotational 
invariance is maintained, thus enabling better comparisons with the continuum theory. 

Even though the main stress of this paper has been on the numerical approxima- 
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tions, the fuzzy disc is also an interesting object from the mathematical point of view, 
and its setting on a more rigorous footing would not only have an intrinsic value, but 
also help understanding the role of the limits in the approximation of field theory and 
renor malizat ion . 



A Conversions between different basis 

In this section we collect some useful conversions formulas between the various basis. 
They are all standard expressions, but care has to be taken because of the nonstandard 
commutation relation (|2.2j) . The resolutions of the identity are 

OO I „ 

1 = £ l"> <*l = - e I \*) (*\ ■ (A- 1 ) 



n=0 



Coherent states are an overcomplete basis and they are not orthogonal, their inner product 
is: 

\z\ 2 + \z'\ 2 -2zz' 

( z \ z') = e-^-e , (A.2) 

while the inner products between a coherent state and an eigenstate of N are given by: 

\z\ n) = e 26 



VnW n 



(n| :) = e_T 7SP ' (A - 3) 



Known is possible to calculate (p mn from 



mi ^ fc) Tay Vime^ 
v>ik = L <Pi- q k- 9 flg -, » ( A - 4 ) 



while the converse is given by 

min(n,m) /_i\p 

P =o p\^(m-p)\(n-p)\O m+n 
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